Ks-band (2.14/im) imaging of southern massive star formation regions 
traced by methanol masers 
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Abstract: 

We present deep, wide-field, Ks-band (2.14 ^m) images towards 87 southern massive star formation 
regions traced by methanol maser emission. Using point-spread function fitting, we generate 2.14^m 
point source catalogues (PSCs) towards each of the regions. For the regions between 10° < I < 350° 
and |6| <1, we match the 2.14 /im sources with the GLIMPSE point source catalogue to generate a 
combined 2.14 to S.Ofim point source catalogue. We provide this data for the astronomical community 
to utilise in studies of the stellar content of embedded clusters. 

Keywords: stars: formation, stars: early type, infrared: stars. Galaxy: stellar content, masers 



1 Introduction 



Stellar populations towards young star formation re- 
gions are deeply embedded within large column den- 
sities of molecular gas and dust. Observations at IR 
wavelengths are particularly powerful for peering into 
these natal cocoons as they suffer much less from ex- 
tinction than optical wavelengths and the young stars 
may also be intrinsically cool. Several excellent large- 
area IR surveys e xist [e.g. 2MASS (|Skrutskie et al.l 
l2006l l. GLIMPSE ('Benja min et al.l 12003')]. the results 
from which have significantly advanced our understand- 
ing of early phases of cluster formation. However, 
these surveys are limited to ~2" resolution which may 
be i nsufficient to resolve t he closest cluster members 
(e.g. iLongmore et al.1 [200^ 1 and the 2MASS Ks-band 
(2.14 /^m) observations are not deep enough to detect 
the most heavily embedded/reddest objects seen in the 
GLIMPSE images. 

In this work we present 2.14/im images and pho- 
tometry for 87 fields associated with massive star for- 
mation regions. These were originally targetted due 
to their I RAS colours in a s e arch for methanol maser 
emission (IWalsh et al.lll997l. [ 199811 and most contain 
hot molecu lar cores ([Longmoire et al.l2007l : lPurcell et al.l 
l2006l . l2009l ). By matching source positions with those 
in existing near-IR (2MASS) and mid-IR (GLIMPSE) 
catalogues, we have created composite 2.14 to 8 ^m 
catalogues which are typically 3-4 magnitudes deeper 
at 2.14/im than 2MASS. Detailed analysis of these re- 
gions requires consideration of completeness and con- 
fusion etc., which is specific to individual fields - a 
result of differing levels of differential extinction and 
source confusion as well as differing sensitivity across 
the fields. However, we calculate a number of metrics 
which describe the quality and quantity of the data 
sets. We provide these catalogues and calibrated im- 
ages as a resource for the astronomical community. 



2 Observations and Data Re- 
duction 

The observations were made using IRISflO (Infrared 
Imager and Spectrometer 2) on the 3.9m Anglo Aus- 
tralian Telescope (AAT) at Siding Spring Observatory 
using the Ks-band (1.982-2.306 /im) filter. 87 regions 
were imaged, as listed in Table [T] The images are 
shown in Figures 7 to 92. The observations were taken 
over several different periods between 2003 and 2006. 
For each source, a 3 x 3 image grid was created with 
~ 1' offsets from the pointing centre. The integra- 
tion time at each of the nine grid positions was around 
one minute, with 9x6s exposures. The initial reduc- 
tion was carried out using a data reduction pipeline 
provided by the Anglo-Australian Observatory. The 
detector bias was removed at readout time using the 
Double Read Mode method. The data were reduced 
using the in-house ORAC-DR pipeline with the 'JIT- 
TER_SELF_FLAT_KEEPBAD' ORAC-DEQ recipe to 
correct for dark current, create flat-field images from 
star free pixels in the 9 jittered source fields and apply 
a bad pixel mask to each image. The pipeline then 
corrects for distortion at the outer edges of the im- 
age caused by the slightly curved focal plane. The 9 
individual fields of 7.7' x 7.7' were finally aligned and 
mosaiced together to give a 9.7' x 9.7' field of view, al- 
though only the inner 5.7' x 5. 7' are covered by all 9 
fields. A typical sensitivity in the inner region of an 
uncrowded field would be expected to be ~19.6mags 
(5a) in good conditions. 



1IRIS2 employs a 1024x1024 Rockwell HAWAII-1 
HgCdTe infrared detector with a plate scale of 0.4486 it 
0.0002" per pixel. 

^see http://www.aao.gov.auJ 

•^see ^http: / /www.oracdr.orgj 
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2.1 Co-ordinate Frame Calibration 

The blind telescope pointing error is estimated to be 
~3" but absolute coordinates were calculated by com- 
paring non-saturated stars to those in the 2MASS cat- 
alogue. Using the koords program in the KARMAQ vi- 
sualisation package, the image coordinates were matched 
to better than 0.1 pixel accuracy. We take the error in 
the absolute coordinates to be the error in the 2MASS 
catalogue of 0.3". 

2.2 Source Extraction 

Sources were extracted from each of the IRIS2 im- 
ages using the daophot tasks in the IRAF0 package 
to produce a 2.14 /im source catalogue (see Section FTll 
for more details). Examination of the residuals af- 
ter removing the stellar flux from the image shows 
the stars have generally been well extracted despite 
different background levels across the image (caused 
by varying extended emission and extinction) and the 
crowded field. Figure [T] shows a histogram of the num- 
ber of sources detected per region (which varied be- 
tween ~6,000 to 25,000, corresponding to stellar den- 
sities ^^60 — 250 stars/arcmin^) and the median source 
photometric fit uncertainty per region ('^O. 06 mag). 

2.3 Photometric Calibration 

The photometry was calibrated using the 2MASS cat- 
alogue towards each region (see Section [3.1(1 . The ex- 
tracted IRIS2 sources (Section 12. 2p were first matched 
with the 2MASS catalogue as outlined in Section [2.41 
The mean offset between the two catalogues was then 
calculated by plotting the difference between the 2MASS 
magnitude and IR1S2 instrumental magnitude of the 
matched stars as a function of the measured IR1S2 in- 
strumental magnitude over a suitable magnitude range 
(typically ^12 to 14, to avoid both 1RIS2 saturation 
and the poorer 2MASS sensitivity). Figure [5] shows a 
histogram of the offset between the IRIS2 instrumen- 
tal magnitude and 2MASS magnitude for each region. 
The range in values is due to the different airmass of 
the IRIS2 observations. The uncertainty in IR1S2 ab- 
solute photometry was estimated from the standard 
deviation between IRIS2 and 2MASS magnitudes over 
a reliable magnitude range (-12"^ to 14*'' magnitude) 
after this offset had been applied. Figure [5] shows a 
histogram of the estimated photometric accuracy to- 
wards each region, which is ~0.13 magnitudes aver- 
aged over all the fields. This is approximately twice 
the statistical error from the photometric fitting (see 
Figure [TJ. Some of this scatter is likely to be due 
to the difference between the MKO photometric sys- 
tem used by IR1S2 and the 2MASS photometric sys- 
terc0 ([Carpenter 2001). Of particular relevance to this 
work, the Ks filters of the two photometric systems 
are responsive over different wavelength ranges - from 
1.982 to 2.306 ^m for IR1S2 and 1.915 to 2.384 ^im 

■'http: / /www.atnf.csiro.au/computing/software/karma/ 
^http:/ /iraf.noao.edu/ 
'^http: / /www. astro. caltech.edu/~jmc/ 
2mass /v3/ transformations / 



for 2MASS. Although converting to the 2MASS pho- 
tometric system may marginally reduce this scatter, 
the colour correction terms are not characterised for 
sources with (J-H) > 1.5, (H-K) > 1 and (J-K) > 
2. These are the reddest, most embedded sources and 
thus likely to be of most interest to studies of star for- 
mation. 

Figure[3]shows histograms of brightest, median and 
faintest magnitude of sources detected per region. De- 
pending on the field, stars brighter than ~11*'' magni- 
tude in the IRIS 2 fields suffer from saturation. Column 
11 in Table[T]gives the minimum reliable magnitude for 
each region, corresponding to the maximum flux per 
field over which the array responded linearly. Photom- 
etry for stars brighter than this should be taken from 
the 2MASS catalogue. 

Table [T] also provides information on the extracted 
IRIS2 2.14/im source populations, image quality and 
photometric calibration towards each of the regions. 
The photometric data themselves are provided in an 
online catalogue in Table [2l We provide here the first 
five lines and a description of the PSC contents. These 
are further combined with the photometry for matched 
sources in the GLIMPSE catalogue (see Section [2.4p . 
The first five lines of the combined catalogue are pro- 
vided here, with the rest in the online material. 

2.4 Matching sources at different wave- 
lengths 

Analysis of source colours relies on the ability to ac- 
curately match sources in the catalogues at different 
wavelengths. With too many sources (>10**) to check 
manually and crowded stellar fields, there is a poten- 
tial problem of either matching a source multiple times 
or mismatching the sources completely. The resulting 
mismatched sources may replicate the large colour dif- 
ferences of embedded sources and hence contaminate 
the sample. To test the matching accuracy, we gen- 
erated and matched synthetic catalogues with known 
offsets from a list of absolute source positions. In or- 
der to ensure a realistic spatial distribution of sources, 
including over densities due to clustering and under 
densities due to extinction, the absolute coordinates 
were taken firom the 2MASS PSC of the region. The 
coordinates recorded in the synthetic catalogues were 
randomly offset from the absolute values by a 2D Gaus- 
sian distribution. Several synthetic catalogues were 
produced using the known absolute pointing error of 
the relevant datasets to generate the positional offsets. 
The sources in the synthetic catalogues were given a 
unique ID and then matched in the same way as the 
observed datasets using the tmatch routine in the ta- 
bles. ttools package in IRAF. The routine matches ev- 
ery source in the first catalogue with the nearest coun- 
terpart in the second catalogue within a user defined 
matching radius. Comparing the ID of the matched 
sources in the synthetic catalogues it is possible to un- 
ambiguously identify which sources have been matched 
correctly. 

The IRIS2 images were all individually registered 
with 2MASS images so have a relative astrometric un- 
certainty of A/fl/s2 = 0.04" (0.1 pixel) [see Section[TT] • 
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Number of sources per region 

Mean = 14571.58 
Standard Deviation = 4779.44 



[IRIS2] median error 
Ivleon = 0.06 
Standard Deviation = 0.02 




10' 1.5x10' 2x10 
Number of sources 



0.05 0.1 0.15 

[iRiS2] medion error (mag) 



Figure 1: [Left] Histogram of the number of sources detected per region and [Right] the median source 
photometric fit uncertainty per region determined using DAOPHOT. 



[IRIS2] - [2IVIASS] offset 
Mean = 0.96 
Standard Deviation = 0.27 



Standard deviation offset from zero point magnitude 
Mean = 0.13 
Standard Deviation = 0.05 




0.5 1 1.5 2 

[IRIS2] - [2MASS] offset (mag) 



0.2 0.4 

Standard deviation offset from zero point magnitude 



Figure 2: Histograms iUustrating photometric cahbration outlined in Section [Left] Histogram of the 
[IRIS2] — [2MASS] instrumental magnitude offset for each region, used to calibrate the IRIS2 catalogues. 
[Right] Histogram of the estimated photometric calibration accuracy towards each region. This was derived 
from the standard deviation between IRIS2 and 2MASS magnitudes of sources over a reliable magnitude 
range (^^12'^ to 14*'^ magnitude) after the calibration factor had been applied. 
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[IRIS2] minimum magnitude [IRIS2] median magnitude [IRIS2] maximum magnitude 

Meon = 8.61 Mean = 15.40 Mean - 18.01 

Standard Deviation = 0.79 Standard Deviation = 0.30 Standard Deviation - 0.64 




[IRIS2] minimum magnitude [IRIS2] median magnitude [IRIS2] maximum magnitude 



Figure 3: Histograms of brightest, median and faintest magnitude of sources detected per region. 



The astrometric uncertainty in the 2MASS and GLIMPSE 
surveys are both quoted as 0.3". However, the GLIMPSE 
survey coordinates were also registered with 2MASS 
data so the relative uncertainty between the IRIS2 
and GLIMPSE catalogues should be considerably less 
than the assumed error of Airis2-glimpse = 0.3". 
Matching two synthetic catalogues generated with an 
uncertainty of A/_r/s2, an equal number of sources and 
a search radius of 2 pixels, produces a 100% recov- 
ery with no confusion due to doubles or mismatches. 
Matching catalogues generated with an uncertainty of 
A/ijjs2 but different numbers of sources also recov- 
ers 100% of the sources correctly but includes ~0.05% 
double sources. With different numbers of sources in 
each catalogue, the matching order then becomes im- 
portant. For example, matching the catalogue with 
larger to smaller numbers of sources may produce 'mul- 
tiple matches', where several sources in the larger cat- 
alogue are within the matching radius of a source in 
the smaller catalogue. However, the number of cor- 
rectly matched sources is the same, irrespective of the 
matching order. Varying the matching radius has lit- 
tle effect as long as it is significantly larger than the 
combined astrometric uncertainty. 

Based on the above analysis, when matching the 
IRIS2 PSC with the GLIMPSE PSC and a matching 
radius of 1.2", we expect to recover at least 99.8% of 
the sources correctly with a maximum of 0.1% doubles 
and 0.12% mismatches. We therefore conclude that 
mismatches should not seriously affect the combined 
IRIS2/GLIMPSE PSCs. 

3 Other Catalogue Data 

3.1 2MASS - 2 micron All Sky Sur- 
vey 

Using two 1.3m telescopes in Arizona and Chile, the 
2MASS survey scanned the entire s ky at the J, H 
and K s bands (1.25, 1.65 and 2.17/im) (|Skrutskie e t al.l 
120061 ;). The survey has a lOcr limiting magnitude of 
15.8, 15.1, 14.3 ± 0.03 mag at J, H and Ks respectively 
with a pixel size of 2.0" and 0.1" pointing accuracy. We 
make use of this data to provide a reference frame for 
the images as well as to provide the flux calibration. 



3.2 GLIMPSE - Galactic Legacy In- 
frared Mid-Plane Survey Extraor- 
dinaire 

Using the infrared array camera (IRAC) on-board the 
Spitzer Space Telescope (SSlQ), the GLIMPSE survey 
has observed the galaxy at 3.6, 4.5, 5.8 and 8.0/im with 
a 1.2" pixel size betwe en 10° < \l\ < 65° and |6| < 1° 
([Benjamin et al.ll2003l ). The photometric accuracy is 
0.2 and 0.3 mag for bands 1/2 and 3/4 respectively. 
The astrometric accuracy of the point source catalogue 
is ~0.3". We combine this with the AAT data to con- 
struct a source catalogue from 2.14 to 8.0/im. 



4 The data 

Figures 7 to 92 show the IRIS2 2.14 /im images to- 
wards each of the regi ons overlayed with the methanol 
maser emission from IWalsh et af] (119981 ). Note that 
in most cases, regions of extinction are readily ap- 
parent. Rarely is there a 2fim point source associ- 
ated with a methanol maser. In many cases, exten- 
sive regions of 2fim nebulosity are seen. Several tar- 
getted regions overlap within a few arcminutes (e.g. 
G14.99-0.70 and G15.03-0.68 in the vicinity of M17) 
but each field was imaged independently to ensure uni- 
form coverage (within the observational systematics) 
in the 5.7' x 5.7' inner region of the mosaic with opti- 
mum sensitivity (see Section [2]|. 

Table [5] shows an example extract from the first 
five lines of the IRIS2 Ks-band (2.14 fim) point source 
catalogue for a sample region. Sections l4.1l fc lT2] exam- 
ine the issues of seeing and completeness in these Ks- 
band (2.14 /im) PSCs including photometry for these 
sources. 

The fits images, IRIS2 PSCs and IRIS2/GLIMPSE 
PSCs for each region are available through the Cen- 
tre de Donnees astronomiques de Strasbourg (CDSfl 
Photometric calibration for the fits images is provided 
by the fits header keyword, "MAGZERO" , which gives 
the magnitude in each image corresponding to a flux 



^http: / /ssc. spitzer. caltech.edu/ 
*http: / /cds. u-strasbg.fr/ 
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count of This is the same value as that given in 
column 6 of Table [T] 



4.1 Seeing 

Figure |4] shows a histogram of the seeing conditions 
during observations of each of the regions. The ma- 
jority of regions were observed in seeing of 1.2 to 1.4 
arcseconds, a noticeable improvement on the ~2" res- 
olution of both 2MASS and GLIMPSE. The other two 
plots in Figure show the effects of seeing on the 
brightest magnitude and number of sources detected. 
When the seeing is poor, the flux from bright stars is 
spread over more pixels and brighter stars can be re- 
covered. Given the crowded nature of the observed 
fields towards the Galactic plane, it is unsurprising 
that the number of sources detected increases dramat- 
ically as the seeing decreases. The inverse linear re- 
lation between source counts and seeing suggests the 
difference in source counts resulting mostly due to the 
seeing present at the time of each observation. 



90 per cent completeness histogram 
Mean = 15.90 
Standard Deviation = 0.£ 




90 per cent completeness 



Figure 5: Histogram of estimated 90% complete- 
ness using a simple, empirical method of calculat- 
ing the turnover in the histogram of the number of 
stars as a function of magnitude, following L07B. 



4.2 Completeness 

iLongmore et al.l (l2007l . hereafter L07B) used artificial 
star recovery to investigate the spatial variation in 
point source sensitivity as a function of wavelength for 
the field associated with G305. 2-1-0.2. Having calcu- 
lated the PSF for each image, a grid of artificial stars 
of the same magnitude separated by 30" was placed 
across the image. The same automated finding tech- 
nique outlined in Section [5] was then used to calculate 
how many of the artificial stars were recovered. By 
shifting the 30"x30" grid of artificial stars in small 
steps through the image, it was possible to measure 
the completeness at 5" intervals without the PSF of in- 
dividual artificial stars overlapping. The process was 
repeated by increasing the artificial star magnitudes 
in steps of 0.5 mag until no more stars were recov- 
ered. By recording the largest recovered magnitude at 
each position and wavelength it was possible to build a 
three-dimensional picture of the point source sensitiv- 
ity across the region. The relative completeness as a 
function of position was then calculated at each wave- 
length by subtracting the median completeness mag- 
nitude at that wavelength from e very position . This 
method is similar to that used bv iGutermuthI (j2005l ) 
in an analysis of completeness limits of Spitzer IRAC 
data. Through analysis of the results from the ar- 
tificial star recovery method, L07B found a simpler, 
empirical method of estimating the 90% completeness 
limit by calculating the turnover in the histogram of 
the number of stars as a function of magnitude which 
yields approximately the same limit as the full analy- 
sis. This metric is provided in Table 1, but we stress 
that a comprehensive analysis of each field would re- 
quire this completeness limit to be determined individ- 
ually for each. Figure [S] presents the histogram of this 
distribution of turnover magnitude. 



^i.e. m = mzp - 2.5 log(f), where f is the number of 
counts determined for a source of interest in the fits images 



5 Applications 

As an example of what the data has the potential 
to provide once a full completeness analysis is con- 
ducted, we show the results towards G305. 2-1-0.2, as 
reported by L07B. Figure|S]shows the slope of the spec- 
tral energy distribution between 2.14 and 4.5 ^m (i.e. 
a2.i4-4.5) of all sources detected towards i) a control 
field using a combined 2MASS/GLIMPSE PSC, u) the 
target field using a combined 2MASS/GLIMPSE PSC 
and ui) the target field using a combined IRIS2/GLIMPSE 
PSC. The spectral index is related to the categories 
of young stellar objects (Class I, II, III), with younger 
sources tending towards larger spectral indexes. In this 
way the PSC provides a measure of the distribution of 
evolutionary states of the embedded cluster seen to- 
wards G305.2-I-0.2. The considerably more sensitive 
2.14 /^m IRIS2 images means much more deeply em- 
bedded (redder) sources are detected with the IRIS2/GLIMPSE 
PSC compared to the 2M ASS/GLIMPSE PSC. This is 
obvious in the bottom panel of Figure [6] as a significant 
increase in the number of sources with a > compared 
to the middle panel. 
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Figure 4: [Left] Histogram of seeing conditions during observations of the regions, [Centre] brightest 
detected magnitude per region vs seeing and [Right] number of sources detected per region vs seeing. 
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7 Appendix 

Table [T] shows the properties of the IRIS2 2.14 ^m 
PSCs for all the observed regions. Tables [2] & [3] show 
the first five lines of a sample IRIS2 2.14 /im PSC and a 
sample combined IRIS2/GLIMPSE PSC, respectively 
with the full tables available on line. Finally, Sec- 
tion [TT] describes in detail how the DAOPHOT pack- 
age (digi.daophot) in TRAP was used to automatically 
find and perform photometry on stars within the IRIS 2 
fields and hence create the point source catalogues. 



www. publish . csiro. au /journals /pasa 



7 



Table 1: Summary of the data characteristics for each of the regions observed with IRIS2. Columns 
one and two give the Galactic and IRAS names of the far-IR sources whose fields have been i maged. 
Column three and four give the R.A. and Dec. of the associated methanol masers from lWalsh et ah l(Ii998). 
Columns five to thirteen show the parameters returned from the automated source extraction algorithm 
using DAOPHOT (see Section mi for more details), azp gives the standard deviation of the difference 
between the measured IRIS2 and 2MASS magnitudes (in the range 10—12 mags). ZP gives the magnitude 
corresponding to 1 flux count in each image. This may be used to calibrate the fits image for each source, 
as described in Section [H Min, Med and Max give the minimum, median and maximum magnitudes of 
sources in each region. Umed gives the median point-source fitting error in each region. Miurez gives the 
minimum reliable IRIS2 magnitude measurement, sources brighter than which are likely to be affected 
by saturation. Columns twelve and thirteen show the 90% completeness magnitude estimated from the 
magnitude histogram and the total number of sources detected in each IRIS2 field, respectively. The final 
column lists the FWHM of the seeing in arcsec. "The dashes in column 2 show sources with no associated 
IRAS point source. ^6305. 20+0.21 has been inclu ded in this table for completeness but the catalogue 
and images can be found in lLongmore et ab l (l2007D . 
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Table 2: Example extract from the first five lines 
of the IRIS2 Ks-band (2.14/im) point source cat- 
alogue towards a sample region. The first two 
columns give the source right ascension and decli- 
nation in decimal degrees. Columns 3 and 4 give 
the magnitude and magnitude uncertainty. The 
full tables will be available on line. 



R.A. 


Dec. 


Mag. 


AMag. 


272.1961 


-20.1837 


16.133 


0.065 


272.2225 


-20.1838 


17.37 


0.144 


272.2219 


-20.1837 


17.108 


0.119 


272.2214 


-20.1837 


17.31 


0.125 


272.2209 


-20.1833 


15.706 


0.041 
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Table 3: Example extract from the first five lines of the online catalogue of the sources matched between 
the IRIS2 Ks-band (2.14 /im) images and the GLIMPSE point source catalogue. The first column gives 
the GLIMPSE ID number. The second and third columns give the source right ascension and declination 
in decimal degrees. Columns 4 and 5 give the magnitude and magnitude error of the sources at Ks 
band. Columns 6 through 13 give the magnitude and magnitude error of the sources at 3.6, 4.5, 5.8 and 
8.0 /im from the GLIMPSE point source catalogue. In this example table, the magnitudes and errors 
have been rounded to 2 decimal places for clarity. A value of 'null' means no source was detected at that 
wavelength. The full table is available online. Note that these PSCs ONLY include sources matched 
in both the GLIMPSE and IRIS2 Ks catalogues. Photometry for sources detected only in the IRIS 
Ks PSC (likely fainter foreground stars) or GLIMPSE PSC (either due to being too faint, or possibly in 
saturated/confused regions) can be extracted from these individual PSCs. Finally, photometry for sources 
in the IRIS2 Ks PSC brighter than the minimum reliable magnitude reported in column 11 of Table [H 
should be extracted from the 2MASS PSC. 



ID 


RA ,12000 


Dec. ,12000 


Ks 


AKs 


3.6 


A3. 6 


4.5 


A4.5 


5 


8 


A5.8 


8 





A8.0 


GLIMPSE 


(degrees) 


(degrees) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(m 


ag) 


(mag) 


(m 


ag) 


(mag) 


SSTGLMCGO 10. 2102-00. 1517 


272.1905 


-20.17984 


12.86 


0.02 


11.90 


0.10 


11.72 


0.09 


11 


79 


0.25 




ill 


null 


SSTGLMC.GO 10. 2213-00. 1721 


272.2152 


-20.18006 


14.02 


0.02 


11.86 


0.07 


11.43 


0.07 


10 


10 


0.13 


10 


83 


0.01 


SSTGLMC.GO 10. 2238-00. 1767 


272.2208 


-20.1801 


12.02 


0.04 


11.51 


0.08 


11.54 


0.09 


11 


35 


0.15 


10 


53 


0.16 


SSTGLMC.GO 10. 21 18-00. 1511 


272.1907 


-20.17816 


13.25 


0.02 


12.04 


0.08 


12.11 


0.13 


11 


71 


0.19 




ill 


null 


SSTGLMC.GO 10. 2431-00. 2087 


272.2606 


-20.17873 


14.06 


0.02 


12.42 


0.08 


12.06 


0.10 




ill 


null 




ill 


null 
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7.1 Creating the point source cata- 
logues 

This section describes how the DAOPHOT package 
(digi.daophot) in IRAF was used to automaticaUy find 
and perform photometry on stars within the IRIS2 
fields and hence create the point source catalogues. 
Description of co-ordinate frame calibration, photo- 
metric calibration and catalogue matching are found 
in Section [21 

Set initial fitting parameters As three seperate 
detector arrays were used over the course of the obser- 
vations care was taken to measure the array specific 
input parameters (eg readout noise, gain, saturation 
value etc.) for each image. The image characteristics 
used by the star finding and photometry tasks were 
calculated using daoedit to pick out 10 bright, non- 
saturated, isolated stars with as little surrounding ex- 
tended emission as possible. By fitting a radial profile 
to the emission, this task then calculates the FWHM, 
SKY (median) and gsky for each of the stars. After 
calculating a sensible average value for each of these 
parameters from stars with good profiles, the values 
are then entered into several 'pars' tasks (which act as 
inputs into the finding/photometry tasks) as follows, 

• epar datapars: 
fwhmpsf = FWHM; 
sigma = crsKY] 

datamin = SKY-(6xasify); 

datamax = "saturation value" (must be less than 

max value ie max 'good' data value); 

• epar centerpars: 

cbox = 2 X FWHM or 5 (whichever is larger) 

• epar fitskypars: 

salgori = mode (for crowded images); 
annulus = 4xFWHM; 
dannulus = 3.5 x FWHM 

• epar photpars: 
aperture = 1 x FWHM 

• epar daopars: 

psfrad = (4xFWHM) + l; 

fitrad = IxFWHM or 3 (whichever is larger) 

Intial star list After calculating the image param- 
eters, the initial (first-pass) star list was created using 
the daofind task to run the search finding algorithm. 
This task was run several times using different detec- 
tion thresholds between 3 and 5cr and compared by eye 
to maximise finding fainter stars and excluding noise. 

Initial photometry Photometry of the stars in 
the first-pass list was calculated using daofind. Stars 
dumped by the task due to poor fits were manually 
checked to ensure no real stars were dropped. In prac- 
tice, the dumped stars were either low signal noise 
spikes or saturated stars. 

Calculating the point spread function The 

task pstselect was used to select 10 bright, non-saturated, 
isolated stars across the image to calculate the PSF. A 



characteristic PSF for the image was then calculated 
from these PSF stars. The robustness of the PSF was 
checked by subtracting the characteristic PSF from the 
10 selected stars and their nearest neighbours and in- 
specting the residuals. In an iterative process, bad 
PSF stars were removed and the procedure repeated 
until good residuals were seen across the image mean- 
ing a realiable characteristic PSF for that image had 
been obtained. 

Extracting the stars and fitting world co- 
ords Once the optimal image parameters and char- 
acteristic PSF had been obtained, the PSF model was 
fit to all the stars using allstar. The residuals after the 
PSF subtraction were carefully checked to ensure both 
bright and faint stars across the image had been suc- 
cessfully subtracted. The star coordinates were con- 
verted from pixels to world coordiates using wcsctran. 



